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ABSTRACT 

We study the emission patterns and light curves of gamma-rays in the pulsar 
magnetosphere with a current-induced magnetic field perturbation. Based on 
the solution of a static dipole with the magnetic field induced by some currents 
(perturbation field) , we derive the solutions of a static as well as a retarded dipole 
with the perturbation field in the Cartesian coordinates. The static (retarded) 
magnetic field can be expressed as the sum of pure static (retarded) dipolar 
magnetic field and the static (retarded) perturbation field. We use the solution 
of the retarded magnetic field to investigate the influence of the perturbation fleld 
on the emission patterns and light curves, and we apply the perturbed solutions 
to calculate the gamma-ray light curves for the case of the Vela pulsar. We find 
out that the perturbation field induced by the currents will change the emission 
patterns and then light curves of gamma-rays, especially for a larger perturbation 
field. Our results indicate that the perturbation field created by the outward- 
fiowing (inward-fiowing) electrons (positrons) can decrease the rotation effect 
on the magnetosphere and makes emission pattern appear to be more smooth 
relative to that of the pure retarded dipole, but the perturbation field created by 
the outward-fiowing (inward-fiowing) positrons (electrons) can make the emission 
pattern less smooth. 

Subject headings: gamma rays: theory - pulsars: general - stars: neutron 
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Introduction 



Since pulsed gamma-ray emission from a large number of pulsars by the Large 
Area Telescope (L AT) on the Fermi Gamma-Ray Space Telescope have been discovered 



flAbdo et al. 



20101), much information are available for constraining pulsar magnetosphere 
geometry. Rotation-powered pulsars are widely believed to have their magnetospheres in 
which charged particles are accelerated to relativistic energy and generate multiband spectra 
of pulsed emission. In order to simulate pulsar magnetosphere, various approximations have 



been proposed. Vacuum magnetosp 



lere model is analytic model for the electromagnetic 



field of a rotating magnetic dipole ( iDeutsch 



19551 ). but it is not an appropriate physical 



model of an active pulsar magnetosphere filled with charges and currents. Force-free 
magnetosphere is probably a closer approximation to a real pulsar than the vacuum solution, 
however it is not a truly self-consistent model since the production of pair plasma requires 
particle acceleration and thus a break down of force-free conditions in some regions of the 
magnetosphere. At present, there are two kinds of three dimensional (3D) emission models: 



one is based on the geometrica 
frame of a retarded dipole (e.g.. 



consideration such as two pole caustic (TPC) model in the 



Dyks fc Rudak 



2003 



Dyks et al. 



2004 



Fang fc Zhandboioh 



and an nular gap (AG) model in the frame of free-force approximation (IBai fc Spitkovsky 



2010a) 



: another is based o n both physical and geometr ical consideration such as slot 



fSG) (iHarding et al. 



Cheng et al. 



2000 



20081) and outer gap mod e- 



Zhang fc Cheng 



2000 



2001 



Romani fc Yadigaroglu. I.-A. 



Zhang fc Cheng 



2002 



Tang et al. 



1995 



2008 



Zhang fc Li 



20091 



Watters et al. 



2009 



Li fc Zhang 



2010 



Wang et al. 



20111 ) 



How we can describe a pulsar magnetosphere filled with charges and currents? Since 
the currents will induc e magnetic field pertu r bation , resulting to a distorted magnetic field 



relative to dipole field. 



Muslimov fc Harding] (120091 ) described an approximate perturbation 



field for a pair-starved current flow in the open zone in the frame of static magnetic dipole. 
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In their treatment, the magnetic field is described as B = B^*^-* + eB*-^\ where B*-'^^ is 
a pure dipole magnetic field anchored into the neutron star (NS), e is the perturbation 
amplitude which determines the strength of the B*-^-* relative to the dipole one , and B^^-* 
is the magnetic field generated by the self-consistent electric currents in the domain of 
the magnetos phere with open field l i nes, w here a singularity exists on the symmetry axis 



(details see in lMuslimov fc HardingI (120091 )). However, although basic features of the pulsar 
magnetosphere in a static dipolar magnetic field approximation can be understood, it is 
more realistic that the magnetosphere be approximated as a rotating inclined magnetic 
dipole. The magnetic field lines of a rotating inclin ed dipole are especially d ifferent from 



those of a static dipole for large inclination angles. 



Romani fc WattersI (120101 ) investigated 



the magnetic field structure with a curren t-induced field in the frame of the retarded 



magne tic dipole. To apply the s olution of 



Muslimov fc HardingI ( l2009l ) to the retarded 



dipole, iRomani fc WattersI ( 120101 ) made following treatments: (1) shifted the magnetic axis 
(r, 6b = 0, 0b) to match that of the swept back dipole (i.e. mapped the magnetic axis line 
onto the swept back curve) and (2) exponentially tapered the singularity appeared on the 
symmetry axis so that the field lines can be integrated to determine the last closed field line 
surface. Obviously, such a treatment method is highly simplified. 

In this paper, we study the emission patterns and light curves in the pulsar 
magnetosphere with a current-induced magnetic field perturbation. We derive the analytic 
solution of the retarded magnet i c field with a current- induced magnetic field; this solution 



is that of 



Muslimov fc Harding 



(120091 ) when rotati ng effects are ignore d . We find that our 
calculated results are different from those given by iRomani fc WattersI (120101 ). for example 
open zone boundary (polar cap) foot-points. Using such a solution, we calculate emission 
patterns and light curves for the Vela pulsar with different inclination angels in the frame of 
two pole outer gap model, where a self-consistent treatment is that both simulation of the 
magnetosphere and relativistic effects are performed in the inertial observer's frame (lOF) 
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[e.g., 



Takata et al. 



2007 



Bai fc Spitkovsky 



2010b|). The paper is organized as follows. In 



§2, the solution of the magnetic field structure with current-induced magnetic field is given. 
Various emission patterns and light curves are calculated in §3. Finally, a brief discussion 
and conclusions are given in §4. 



2. Magnetic Field Structure 



In a pulsar magnetosphere, radiating charges accelerated in gaps will create some 
currents in the open zone. If gap-closing pair front produce densities comparable to the 
co-rotatio n value for at least sorne of the open field lines, this cu rrent will induce a magnetic 



field (e.g. 



Muslimov &: Harding 



2009 



Romani &: Watters 



local magnetic field can be described as ( jMuslimov fc Harding 



20101). Fo r a static dipole, the 



20091) 



B = B(^) + eB(^) , 

where B^'^-' is the pure dipole magnetic field and is given by 

1 



B(d) 



[3( 



m ■ r r — m 



(2) 



where m is the magnetic moment and its three components in Cartesian coordinates can 
be expressed as (the derivation is shown in Appendix A) 



my = - [3yx5f + (3i/2-2r2)5f ] 
m, = - [2xzB^^^ + 2yzBf^ + (3^' - 2r^)B^^^] 



(3) 



where r = a/x^ + + z^, and B^\ b'^\ and B^f^ are three components of the magnetic 
field given by Eq. ([2]). In Eq. (111), e is the perturbation amplitude, and B*-^^ is the magnetic 
field generated by the self-consistent electric currents in the domain of the magnetosphere 
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with open field lines. After introducing cosine of the angle between the neutron star(NS) 
rotation axis and radius- vector of a given poin t g = cos a cos 6^ + sin a sin 6' cos and its 



derivatives over 9 and 



Muslimov fc Harding (120091 ) derived the analytic expressions of 



Eq. ([T]) in magnetic polar coordinates for the static dipole, where s = 1 determines the 
symmetry axis in magnetic coordinates which is the NS rotation axis. In Appendix A, we 
change expressions (5^, Bq, B^) of Eq. ([T]) in magnetic polar coordinates to those {B^^ By, 
Bz) in Cartesian coordinates, and then give three components of the magnetic moment m' 
in which the magnetic field is given by Eq. ([1]) (i.e. the perturbation field is included) by 
changing {B^\ By^\ Bi^^) to (5^., By, B^) in Eq. (3). In such a treatment, the magnetic 
moment m in the pure dipole magnetic field changes to m' in the magnetic field given by 
Eq. ([1]). We call m' as an effective magnetic moment. The difference between m and m' 
depends on the perturbation amplitude e, (5m = |m — m'| =0 when e = 0. 

Since the perturbatio n field becomes singular along the rotation axis where s — )■ ±1, 



Romani &: WattersI (120101 ) exponentially tapered this singularity (oc 1 — exp[(|s| — \)/as], 
with as = 0.25) so that the field lines can be integrated to determine the last closed field 
line surface. Here we will use this method. 

The structure of the magnetosphere can still only be described by a specific features 
such as the last closed field lines. The footprints of the last closed field lines on the stellar 
surface define the polar cap shape. In the static approximation, the polar cap shape 
is a circle with radius of Rpc = Rq{Rq/ Ri.Y^'^ for an aligned dipole, where Ri^ = c^l^^ 
is the light cylinder radius, and Rq is radius of the neutron star, we define the initial 
polar cap's edge using {xo,yo,zo) = [-RpcCos(0p), i?pcsin(0p), - i^^J ], where 0p is 
the azimuthal angle about the magnetic axis. Then we try to find the scaling factors qq 
such that (x'n, y'n. z'n) = [aoXo,aoyo, (R^ — alR^y^"^ ] by using the Runge-Kutta method 



(IZhang &: Cheng 



200ll ). Next we define another scaling factor ai to denote the footprints of 
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Fig. 1. — Changes of open zone boundary (polar cap) foot-points with phase for the a = 70° 
static dipole with current-induced perturbations e = 0, ±0.05, and ±0.1. The Vela pulsar's 
parameters are used. 
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Fig. 2. — Changes of open zone boundary (polar cap) foot-points with phase for the a = 70° 
retarded dipole with current-induced perturbations e = 0, ±0.05, and ±0.1. The Vela pulsar's 
parameters are used. 
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Fig. 3. — Changes of open zone boundary (polar cap) foot-points with phase for the a = 70° 
retarded dipole with e = for the periods of Crab, Vela, and Geminga pulsars. 
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Fig. 4. — A three-dimensional view of last closed field lines for a retarded dipole with or 
without the perturbation field inclined to the rotation axis with angle a = 50° for Vela 
pulsar. The magnetic field is given by Eq. (|5]). The panels from top to bottom represent 
that the cases of e = 0.1, e = 0.0, and e = —0.1. For each panel, the thick red line marked by 
A represents the magnetic axis of pure static dipole, and the thick green curved line marked 
by B represents the magnetic axis of the retarded magnetic field given by Eq.Q. 
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the open field lines as (x, z) = {uiXq, oit/q, [Rl — (x^ + y^)]^/^). As an example, we show 
the variation of open zone boundary (polar cap) foot-points with phase for the a = 70° 
static dipole with current- induced perturbations e = 0, ±0.05, and ±0.1 for Vela pulsar 
parameters (i.e. its period is P = 0.0893 s and surface magnetic field is B = 3.3 x 10^^ G) 
in Fig. [H It can be seen that the open zone boundary appears two symmetric peaks if the 
current- induced B perturbation is not included (i.e. e = 0). However, such a symmetry 
breaks when the current-induced B perturbations is included (i.e. e 7^ 0). 

For an inclined dipole of constant magnitude, rotating with angular frequency Q about 
the z-axis in which the magnetic field consists of pure retarded dipolar magnetic field and 
perturbation field, the magnetic moment is 



m = m'^ ± em 



p ' 



(4) 



where m'' and nip are the magnetic moments of pure retarded dipolar magnetic field B*" 
and perturbation magnetic field Bp, respectively and the expression are given in Appendix 
B. In such a case, the magnetic field is 



B 



B"^ ± eB: 



(5) 



where the expressions of B"" and Bp are given in Appendix B. 

In Fig. m we show the variation of open zone boundary foot-points with phase for 
the a = 70° retarded dipole with current-induced per t urbat ions e = 0, ±0.05, and ±0.1 



for comparison with the results of 



R,omani fc WattersI (I2OI0I ). Since the rotating effect in 



the retarded dipole is included, two symmetric peaks disappear even the current-induced 
B pert urbation is not included (i .e. e = 0). This result with e = is consistent with 



that of iRomani fc WattersI (I2OIOI ) (see their to p panel of Fig. 2), but th ere is a small 



difference between our result and the result of 



Romani fc WattersI fl2010f ). In fact, the 



shape of the open zone boundary foot-points depends on not only the magnetic inclination 
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angle but also the pulsar's period, we show the shapes for the periods of Crab, Vela, and 
Geminga with a = 70° in Fig. E] to account for this fact. Mo reover, the cylindrical radius 
r± = 1.2i?L = 1.2c/n was used in iRomani &: Watterd (120101 ). but r± = Ri^ is used in 



our treatment, leading to the difference. With increasing B perturbation (i.e. increasing 
e), open zone boundary is significantly different from that with e = 0, as shown in Fig 



El Mo re importantly, for the singulari ty in the pair-starved field of 



Muslimov fc Harding 



(120091 ) . we use the smooth method of iRomani &: Watterd (l2010l ) to deal with it in static 
dipole field (see Fig . [Tl), s o this singularity does not exist in retarded dipole field. However, 



Romani &: Watterd (120101 ) took the solution of 



Muslimov &: Harding] (120091 ). and shifted 



the magnetic axis to match the retarded dipole by simply mapped the magnetic axis line 



onto the swept back curve. There 
given by 



"ore, our results with e 7^ are very different from those 
Romani fc Watterd (j2010l ). Physically, the perturbation fields with e < and e > 
should have opposite roles on pure retarded dipole field, our results show this property, so 
we believe that our results are reasonable. 

In the retarded dipole, the magnetic axis is curved relative to that of pure static dipole. 
In Fig. m we show a three-dimensional view of last closed field lines for a retarded dipole 
with or without the perturbation field inclined to the rotation axis with angle a = 50°, 
where the thick red line marked by A represents the magnetic axis of pure static dipole, 
and the thick green curved line marked by B represents the magnetic axis of the retarded 
dipole with the perturbation field (e 7^ 0). It can be seen that the perturbation field has an 
important role in both the shape of the last closed field lines and the magnetic axis. For 
example, when e > 0, the shape of the last closed field lines becomes more smooth and the 
magnetic axis is more curved on the side of the magnetic axis of pure retarded dipole (see 
top panel of Fig. H]). 
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3. Emission Patterns and Light Curves 

In order to determine photon emission from the outer gaps, 3D pulsar's magnetosphere 
should be simulated and relativistic effects including photon abberation and time-of-fiight 
phase shifts need to be taken into account. The geometry of the photon emission is usually 
expressed as {(, $), where ( is the polar angle from the rotation axis and $ is the phase 
of rotation of the star. Here we consider the photon emissions from two poles of pulsar 
magnetosphere and the photon emission geometries in the inertial observer's frame (lOF). 



The particle motion in the lOF is described by (ITakata et al. 



20071 ) 



n = [l3o cos (p{r)]h + [/3osinv9(r)]b_L + /ScoG^, (6) 

where ip{r) is the pitch angle. The third term represent corotation with the star, Pco = p^/c. 
The quantity /3o at each point is determined by the condition that |n| = 1. The unit vector 
h± perpendicular to the magnetic field line is 

= ±[{cos5(f))k + {smS(j))k x b], (7) 

where ± corresponds to the gyration of the positrons (+) or electrons (-), 6(j) refers to 
the phase of the gyration, and k = (b ■ V)b/|(b • V)b| represents the unit vector of the 
curvature of the magnetic field line. Therefore, the emission direction (C, $) in the lOF is 

cosC = n^, (8) 

and 

$ = -$n-r-n, (9) 

where — $„ is the azimuthal angle of the emission direction, and r is the emitting location 
in units of the light radius. 

We now consider the emission patterns and the examples are shown in Fig. [51 In this 
figure, the emission patterns with e = 0.1, 0, and -0.1 for a = 50°, 60°, and 70° are given. 
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Fig. 5. — Emission projections onto the (C, $)-plane for various magnetic inclination angles 
and for e = 0.1, 0, and -0.1. The black lines represent the emission from one pole and the 
grey lines represent the emission from another pole. Vela pulsar parameters are used. 
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Fig. 6. — Radial distances to the null charge surface for the last open field lines for a = 50' 
and different values of e. The Vela pulsar's parameters are used. 
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Fig. 7. — Light Curves for different magnetic inclinations, view angles, and e. Panel A: 
a = 50°, C = 60°, 75°, 90°, and e = -0.1, 0, 0.1. Panel B: a = 60°, C = 60°, 75°, 90°, 
and e = -0.1, 0, 0.1. Panel C: a = 70°, C = 60°, 75°, 90°, and e = -0.1, 0, 0.1. Both 
outward emission from one pole and inward emission from another pole are included and the 
azimuthal range of outer gaps is hmited by using the method of Fang & Zhang (2010). The 
Vela pulsar's parameters are used. 
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For the magnetic field with the perturbation field, the emission patterns is different from 
those for the pure retarded dipole. The emission pattern for e > appears to be more 
smooth relative to that of the pure retarded dipole, but the emission pattern for e < is 
less smooth (see Fig. |5]). Note that the differences depend on the value of e. 

Before calculating the light curves, we conside r the change of t he nu ll charge surface 



with the phase. In the two-pole caustic model of 



Dyks &: RudakI ( 120031 ) . photons are 



emitted uniformly in the gap, and the gap extending 



is confined to the last open field lines. 



Fang fc Zhang 



rom t 



le star surface to high altitudes 



( I2OIOI ) made a revised version of the 



two-pole caustic model, they pointed out that although acceleration gaps can extend from 
the star surface to the light cylinder along near the last open field lines, the extension of 
the gaps along the azimuthal direction is limited because of photon-photon pair production 
process. In such gaps, high-energy photons are emitted uniformly and tangentially to the 
field lines but cannot be efficiently produced along these field lines where the distances to 
the null charge surface are larger than ~ 0.9 times of the distance of the light cylinder. 
In Fig. El we show the changes of the null charge surface with the phase for a = 50° 
but different values of e. It can be seen from this figure that there are substantially 
difference between the null charge surfaces of the pure retarded dipole with and without 
the perturbation field, and the difference becomes larger when e increases. 



Fang fc ZhaneJ f l2010h 



In calculating the light curves, we use the method given by 
i.e. the emission region is limited in the range of TjiuuZ-Rl < 0.9 and Pmax/-RL = 0.95, 
where P rna^ is the distan c e from the rotation axis. On the other hand, we use the method 



given by iDyks fc RudakI (j2003l ) to smooth the light curves, i.e. the model light curve was 
calculated for electrons distributed evenly along the polar cap rim; the density profile across 
the rim was assumed to be the Gaussian function /(^m) symmetrical about 6^ = 6^^, with 
a = 0.0256'pc, where 6^ is the magnetic colatitude of magnetic field lines' footprints at the 
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Fig. 8. — Comparison of light curves (black solid lines) predicted in two pole outer gap 
model with observed light c urve (red lines wit h error bars) in energy region from 0.1 GeV 



to 10 GeV of Vela pulsar (lAbdo et al. 



20091 ). The predicted light curves in upper and 



bottom panels are calculated in the retarded field without any perturbation and with a 
perturbation of e = —0.1, respectively. Note that the predicted light curves come mainly 
from the contribution of one pole since the contribution (grey lines) of another pole can be 
neglected when the values of a and ( listed in the figure are used. 
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star surface, and 9pc ~ ('"ns/''"ic)^^^ is the magnetic colatitude of the rim. We show the 



hght curves in different magnetic inchnations and view angles for e = —0.1, 0.0, and 0.1 in 
Fig. [71 where we have assumed that the inner boundary of the outer gap is located on the 
stellar surface. It can be seen that the change of the light curves with e 7^ is remarkable 
relative to those with e = 0, particularly for larger value of e. Compared to the light curves 
with e = 0, the light curves with e 7^ have following features: (1) peak shapes is changed 
although the peak separations are roughly the same, and (2) the light curves with e > are 
more smoothing but those with e < become more complicated. In Fig. |8l we show the 
light curves predicted in the frame of two pole outer gap model and c ompare them with 



2OO9I). 



observed one of energy region from 0.1 GeV to 10 GeV of Vela pulsar (lAbdo et al. 
The predicted light curve in the upper panel is calculated in the retarded field without any 
perturbation, while the lower panel is calculated in the retarded field with a perturbation of 
e = —0.1. In our calculations for FiglHl the inner boundary of the outer gap is assumed to 
be the null charge surface, the predicted light curves are mainly produced by the emission 
of one pole although the contributions of two poles are included when the values of a anc 



C liste d in the figure are used. Note that the light curve s predicted in 



(I2OIOI ) are based on their one pole outer model given by 



Romani 



Romani fc Watters 



( Il996l ). we can reproduce 



their result with e = in this one pole outer gap model when we used the parameters (i.e. 



a 



72° and ^ = 64°) listed in Fig. 5 of 



Romani &: WattersI (l2010[ l. Comparing our results 



with their results we reached the same conclusion, i.e. that the light curve with e < is 
more consistent with the observed one for the Vela pulsar. 



4. Discussion and Conclusions 



Since some currents in the open zone will be created by the r adiating charges acceleratec 
in the gaps and the currents will induce the magnetic field (e.g., iMuslimov fc Hardingll2009l : 
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Romani &: Wattera 120101 ). it is important to study the structure of the magnetic field with 
the perturbation field, emission patterns, and light curves in pulsar magnetosphere. In this 
paper, we derive the solution of the static (retarded) magnetic field with the perturbation 
field in the Cartesian coordinates, which can be expressed as the sum of pure static 
(retarded) dipolar magnetic field and the perturbation field in in the Cartesian coordinates. 
We have confirmed the reliability of our results by using the solution of the static magnetic 
field with the perturbation field (see Fig. [1]), and then we have used the solution of the 
retarded magnetic field with the perturbation field to investigate the emission patterns and 
light curves in the pulsar magnetosphere (the parameters of the Vela pulsar are used). Our 
results show that the photon emission pattern and light curves for the retarded magnetic 
field with the perturbation field are changed with respective to those for pure retarded 
dipolar magnetic field (see Figs. |S]and[7]). 



Romani Sz Watterd (120101 ) have studied the emission patterns and light curves of the 
retarded magnetic field with the perturbation field by using a highly simplified method. 
However, our method is totally different from their method and the results obtained for 
these two methods are different (see Figs. [2] and [8]). It should be n oted that a more 



realist ic field to describe pulsar magnetosphere is the force free field (IBai fc Spitkovsky 



2010bl ). but such a filled magnetosphere lac ks the acceleration fields required to produce 



powerful 7-rays. 



Romani fc Watterd (120101 ) pointed out that the force free field can be 



approximated as a vacuum dipole field possibly with current-induced perturbation field. 



Fermi are needed in the frame of m ore realistic models such as slot gap (e.g., 



2008; 


Harding & A/ 


uslimov 


2011 


) and outer 


1998; 


Zhang et al. 


2004; 


Hirotani 


2008 


) mo 



Romani 



1996 



(e.g.. 


Harding et al. 


Zhane: & CheneJ 



20081 ) models in which both physical and geometrical 



consideration are taken into account. 
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Appendix 



A. The Static Dipole Field 



Three components of the static dipole field [Eq. ([2])] in Cartesian coordinates are 



' = —{x nix + xyniy + zxriiz) j- 

B^y'^ = ^{yxm^ + y^niy + yzm^) - ^ 



y 

Bi"^ 



-(zxrUx + zyniy + z'^niz) — 



(Al) 



From this set of equations, we can derive three components {rrix, my, ruz) of magnetic 
moment m which is given in Eq. (3). After including the perturbation field B'^-'-^ the 
magnetic moment m i s cha nged to m' = m + enip and the form of Eq. (3) is still valid. 
Muslimov fc HardingI (120091 ) derived the analytic expressions of Eq. ([1]) in magnetic polar 
coordinates for the static dipole, which are 



Bj. 
Be 

Ba, 



B. 



(d) 



COS 9 + e X 



B, 



(d) 



2 ^^"^+1-.^ 



1 - s\ 
ex f ds 



1 



1 ds 



09 ^/T^ sin 9 



(A2) 



B, 



(d) 



e X 



1 ds 



1 



ds\ 



r/3 i_s2 ysm9 d(l) ^^T^ 09 J ' 

where rj = r/R, R is NS radius; x = ^/Rh and i?L is the light cylinder radius; 

s = cos a cos 9 + sin a sin 9 cos </> and a is the angle between the rotation axis and magnetic 

axis. Since the magnetic field in Eq. (A2) are expressed in magnetic polar coordinates, the 
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transformation between magnetic polar coordinates and Cartesian coordinates is 
Br — sin^cos0 Xm + sin^sin0 ym + COS^ Zm , 

ee — sin(6' + 7r/2)cos0 Xm + sin(6' + 7r/2)sin0 + cos(6' + 7r/2) Zm , (A3) 
= cos(0 + 7r/2) x^ + sin(0 + 7r/2) , 

where Zm is parallel to the magnetic axis. Therefore, the components of the magnetic field 
in the Cartesian coordinates are given by 

B^^ = sin e cos (l)Br + sm{e + ■k/2) cos + sin(6' + 7r/2)B^ , 

By^ = sin e sin + sin (61 + 7r/2) sin + sin(0 + 7r/2)B^ , (A4) 

B,^ = cos eBr + cos(^ + t:/2)B0 . 

In principle, converting above components of the magnetic field in the magnetic polar 
coordinates into those in the Cartesian coordinates where the spin axis of the dipole is 
assumed to be along the z-axis, we can obtain the expressions of the magnetic field (Eq. 
(1)) in the Cartesian coordinates. However, since the perturbation field B^-^^ will make 
the magnetic axis to curve, we use (cto, 0o) to describe the deflection of the magnetic 
axis in both magnetic inclination and azimuthal directions. Through fitting the results of 
numerical simulation, we find that ckq and 0o are approximated as 

Xi = -0.02391 - 0.07578Q; + O.lOSTSa^ - 0.03825a;^ , 
X2 = 0.00189 - 0.08582e-2-^^^^^" , 

where r^y is the distance to the spin axis. After eliminating the influence of curved magnetic 
axis, we can use following expressions to obtain the expressions of the magnetic field ( Eq. 
(1)) in the Cartesian coordinates: 

B^ = cos(q;)S^^ + sin(Q;)S^^ , 
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By = By^, (A6) 

Replacing By'^\ and i?!'^'' into -B^;, By, and Bz in Eq. (3), three components of magnetic 
moment m' of the magnetic field given by Eq. (1) can be obtained as follow: 

m'^ = - [{3x^ - 2r^)B^ + 3xyBy + SxzB^] , 

m; = ^[3yxB, + {3y''-2r'')By + yzB,] , (A7) 
K = t; [ixzB^ + 3yzBy + {3z^ - 2r^)Bz] , 

where the perturbation magnetic moment nip can be expressed as 

m^x = - [(3x^ - 2r^)S(^) + 3xyBf^ + 3xzB^}^\ , 

= ^ [3^x4^) + (3y^ - 2r2)sW + t/^Sf^)] , (A8) 
m^, = ^[3xzi?W+3yzS« + (3z2-2r2)5«] . 



B. The Retarded Dipole Field 

B.l. Retarded Magnetic Moment 

For an inclined dipole without the perturbation field of constant magnitude, rotating 
with angular frequency O about the 5-axis, the magnetic moment is expressed as 

m"" = m(sinQ;cos(Ot) x + sinQ;sin(Ot) y + cosa z) , (Bl) 

i.e. the magnetic moment for the static dipole rotates along the rotational direction by 

= (ro — r)/RL- For the same reason, three components for rotating perturbation 
magnetic moment nip are given by 

rrf — rupx cos(Qf) + rupy sin(Qf) , 
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rripy = rupx sm{Qt) + nipy cos (Qt) , (B2) 

= rrip^ . (B3) 

Therefore, the magnetic moment which includes the perturbation field is given by 

m''" = m'" + enip . (B4) 



B.2. Magnetic field with the Perturbation Field 



For the pure retarded dipole, the magnetic field is expressed as (ICheng et al. 



20001 ) 



r ■ 



+ 



+ 

cr^ c^r 



^3 



(B5) 



where r is radial distance, r is the radial unit vector, and c is the light speed. Replacing m 
into m'"" = m"" + enip in above equation, the magnetic field with the perturbation field is 
given by 

B"^ = B'' + eB 



p ' 



where the rotating perturbation field Bp is given by 



Snip 3mp nip 



cr^ 



c^r 



m" m' 

,3 ^2j. 



The three components in the Cartesian coordinates can expressed as 



R — r R — r 
) + Dxsm(^ — ) + SxzTJipz 



Bpx 


= X 


Bp = 


1 


CxCos( 




= y 


Bp = 


1 


CyCOS( 




= z ■ 


Bp = 


1 


CzCos(- 



Rl 
R-r 

Rl ' 



R-r 

+ Dysin(— — ) + Syzrupz 
tir, 



R — r 

+ D,sin(— — ) + (r^ - ?,z^)m 
Rt, 



pz 



where 



(B6) 



(B7) 



(B8) 



Cx rflpx XyUlpy X Tflpx) 2 



ir'^m. 



py - 'iXyVflpx 3X Tflpy^ 
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T 

Dx = {r^vupy - xynipx - x^mpy)—^ + {r^rupx - 'ixyrupy - 3x'^mpx)— 

(r Tfipy Sxyuipx 3x uipy^, 

Cy = {r'^nipy - xyrupx - y'^rUpy)—^ - [r'^rUpx - SxyrUpy - Sy'^rUpx)- 

{r'^rupy - Sxyrupx - ^y'^nipy), 

Dy = {r'^TUpx - xyrupy - y'^rupx)-^ + {r^rnpy - SxyrUp^ - Sy'^mpy)— 



{r'^rupx - ^xyrUpy - Sy'^nipx), 



^2 ^ 



Cz = {-xzuipx - yzrupy)—^ + {Sxzrupy + Syznipx) — + {SxzrUpx + SyzrUpy), 



,^ , ^ , 

^2 ^ 



Dz = (-xzuipy - yzrupx)-^ - (Sxznipx + Syznipy)— + (Sxznipy + SyzrUpx). 



-26- 

REFERENCES 

Abdo, A. A. et al. 2009, ApJ, 696, 1084 

Abdo, A. A. et al.2010, ApJS, 187, 460 

Bai, X. -N. & Spitkovsky, A. 2010a, ApJ, 715, 1282 

Bai, X. -N. & Spitkovsky, A. 2010b, ApJ, 715, 1270 

Cheng, K. S., Ruderman, M. A., & Zhang, L. 2000 ApJ, 537, 964 

Deutsch, Arnim J. 1955, AnAp, 18, ID 

Dyks, J., & Rudak, B. 2003, ApJ, 598, 1201 

Dyks, J., Harding, A. K., Rudak, B. 2004, ApJ, 606, 1124 

Fang, J., & Zhang, L. 2010, ApJ, 709, 605 

Harding, A. K., Stern, J. V., Dyks, J., Prackowiak, M. 2008, ApJ, 680, 1378 

Harding, A. K., Mushmov, A. G. 2011, ApJ, 726, LIO 

Hirotani, K. 2008, ApJ, 688, L25 

Li, X. & Zhang, L. 2010, ApJ, 725, 2225 

Mushmov, A. G. & Harding, A. K. 2009, ApJ, 692, 140 

Romani, R. W. & Yadigaroglu, I.-A., 1995, ApJ, 438, 314 

Romani, R. W. 1996, ApJ, 470, 469 

Romani, R. W., & Watters K. R 2010, ApJ, 714, 810 

Takata, J., Chang, H.-K., & Cheng, K. S. 2007, ApJ, 656, 1044 



Wang, Y., Takata, J., Chang, H.-K., & Cheng, K. S. 2011, MNRAS, tmp 575, 
laFXw: 1102.44741 

Tang, A. P. S., Takata, T., Jia, J. J., & Cheng, K. S. 2008, ApJ, 676, 562 
Watters K. P., Romani, R. W. Weltevrede, P. & Johnston, S. 2009 ApJ, 695, 1289 
Zhang, L., & Cheng, K. S. 1998, MNRAS, 294, 177 
Zhang, L., & Cheng, K. S. 2000, A&A, 363, 575 
Zhang, L., & Cheng, K. S. 2001, MNRAS, 320, 477 
Zhang, L. & Cheng, K. S. 2002, ApJ, 569, 872 

Zhang, L., Cheng, K. S., Jiang, Z. J., & Leung, P. 2004, ApJ, 604, 317 
Zhang, L., & Li, X. 2009, ApJ, 707, L169 



This manuscript was prepared with the A AS L^TgX macros v5.2. 



